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Abstract 

In preparation for an experimental study of magnetorotational instability 
(MRI) in liquid metal, we present non-ideal axisymmetric magnetohydrodynamic 
simulations of the nonlinear evolution of MRI in the experimental geometry. The 
simulations adopt fully insulating boundary conditions. No-slip conditions are 
imposed at all boundaries. A clear linear phase is observed with reduced linear 
growth rate. MRI results in an inflowing "jet" near the midplane and enhances 
the angular momentum transport at saturation. 



Subject headings: accretion, accretion disks 
cal — MHD 



instabilities — methods: numeri- 



1. Introduction 

The magnetorotational instability (MRI) is pro bably the main cause of turbulence and 
accretion in sufficiently ionized astrophysical disks (IBalbus fc Hawleylll998l ). which has in- 
spired searches for MRI in Taylor- Couette flo w. Experiments on magnetized Coue tte flow 
aiming to observe MRI have been performed (ISisan et al.l 12004 ; IStefani et al.l 120061 ) . Those 
experiments have demonstrat ed MRI-like mode s, but not yet for background flows that 
approximate Keplerian disks. ISisan et al.l (120041 ) observed nonaxisymmetric modes and en- 
hanced transport i n a spherical e xperiment, but with a hydrodynamically li n early unstable 
background state. IStefani et al.l (120061 ); iRiidiger et al.l (120061 ); IStefani et al.l (120071 ) demon- 
strated unsteady behavior in a very low-Rm experiment with a strong toroidal as well as axial 
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field; the relationship between the experimental results and linear instabilities is still con- 



troversial (Liu et a 



Szklarski 2007 



2006aJ; iPriede et all 120071 : iLiu et all 120071 iRiidiger fc Hollerbachl 12007 



Liull2008al ). but if these are indeed the linear modes, then the correspond- 



ing ratio of toroidal to axial field in a thin (thickness Ah <C r) disk would have to be 
~ r/Ah, where r is the radius of the a ccretion disk. Some other experiments have been 



proposed or are still under construction (INoguchi et al.l 120021 ; IVelikhov et al.l 120061 ; I Ji et al. 
200ll ; iGoodman fc JillioO^ I. 



Standard MRI modes will not grow unless both the rotation period and the Alfven 
crossing time are shorter than the timescale for magnetic diffusion. This requires that both 
the magnetic Reynolds number Re m = r2 1 r 1 (r 2 — ri)/rj Ga _ (see definitions of rj Ga _, r x and 
r 2 in Fig. [T]and their values in Table. []JJ and the Lundquist number S = V^ 2 (r 2 — ri)/r] Gsi be 
> 1, where V% z = -Boo/ V^PGa is the Alfven speed, pea is the density of liquid Gallium and 

is the background magnetic field parallel to the angular velocity. Because the magnetic 
Prandtl number Pr m = u/r] ~ 10~ 5 — 10~ 6 in liquid metals, where v G& is the kinematic 
viscosity of liquid gallium, Re > 10 6 and fields of several kilogauss for standard MRI must 
be achieved in typical experimental geometries. Recent linear analyses have shown that 
MRI-like mode could grow with much reduced magnetic Reyn olds number and Lundquist 



number in the pres ence of a helical background magnetic field ( IHollerbach &: Riidigerl 12005 



Riidiger et al.ll2005l ). at least for cylinders of infinite or periodic axial extent (Liu et al 2006b). 
The simulations of this paper, however, are limited to axial background fields 

The experiment is complicate d by the large (> 10 6 ) Reyn olds number and by Ekman 
circulation and Stewartson layers (IHollerbach fc Fournierl 12004 ). even though the Princeton 
MRI experimental apparatus has been co nstructed to minimize the circulation by the use 



of ind ependently controlled split endcaps (IKageyama et al.ll2004j ; iBurin et al.ll2006l ; I Ji et al. 



20061 ) . It is known that Ekman circulation is significantly modified when the Elsasser number 
A exceeds unity: A = B^ / (87rp G a?7Ga^) ^ 1 where = v^i^2 is the characteristic rotation 
frequency (jGilmanlll97l[ h At 100% of the maximum designed rotation rate of our experiment 
(TablefjjJ and with I v = 1000 A, where 1^ is the external coil current (Fig. [1]), the Elsasser 
number A = 0.36. 

There have been few published studies o f the nonlinear satu r ation o f MRI in an exper - 
imental Taylor- Couette geometry apart from iKnobloch fc Juried (120051 ); ILiu et al.l (j2006bl ); 
Umurhan et al.l (2007b. af), and even fewer with realistic experimental boundary conditions. 



As noted in ILiu et al.l (j2006bl ). the nature of saturation of MRI in a Couette flow is es- 
sentially diffe rent from accretion disks, in which MRI is b elieved to saturate by turbulent 
reconnection (IFleming et al.l |2000| ; ISano &: Inutsukal 1200 ll ). In accretion disks, differential 



rotation arises from a balance between the gravitational attraction of the accreting body 
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and centrifugal force. Thermal and turbulent energies are probably small compared to or- 
bital ones, so saturation is not likely to occur by modification of the mean flow profile. In 
experiments, however, differential rotation is imposed by viscous or other weak forces, so 
that mildly nonlinear MRI might well change the mean rotation profile. 



This paper is the second of a series, following iLiu et al.l (l2006bl ). Several idealizations 
made in that paper, notably the vertically periodic boundary conditions, are dispensed with 
here. We discuss the linear phase and nonlinear saturation of MRI in finite cylinders with re- 
alistic fluid and magnetic bou ndary conditions. A ll simulations reported here were performed 
with the ZEUS-MP 2.0 code (jHayes et al.l 120061 ) . which is a time-explicit, compressible, as- 
trophysical ideal MHD parallel 3D code, to which we have added viscosity, resistivity (with 
subcycling to reduce the c ost of the induct ion equation) for axisymmetric flows in cylin- 
drical coordinates (r,<p,z) (ILiu et al.l l2006bl ). The imple mentation of full y insulating and 
partially conducting boundary conditions are discussed in ILiu et al.l (120071 ). The computa- 
tional domain is shown in Fig. [1] and the parameters are summarized in Table [H Six coils 
with dimensions as shown were used, with 67 turns in the two coils nearest the midplane 
and 72 in the rest. They are split into two sets of three in parallel, with the upper three 
in series and the bottom three also in series. Currents 1^ were adjusted according to the 
experimental values. Note that in the simulations the magnetic diffusivity r/Ga is fixed to 
the experimental value rjQ a = 2,430 cm 2 s _1 (TableJT]), however the kinematic viscosity is 
varied for the purpose of extrapolation from numerically tractable to experimentally realistic 
Reynolds numbers. 

The measured current waveform is displayed in Fig. [2] for target currents of 1000 A and 
400 A. The waveform displayed has an overshoot at the early stage, a linear decline and then 
a linear ramp where the controller tries to adjust the output voltage to reach the programmed 
set point. This behavior is peculiar to the high current runs (left panel of Fig. [2]). Lower 
currents have a much flatter waveform (right panel of Fig. [2]). In the simulation a wave form 
like that in Fig. [3] is used to approximate the experimental coil currents with ramp time 
t C oii = 0.2 s. 

Since the container is made of stainless steel, not a perfectly insulator, the radial mag- 
netic angular momentum flux {—rB tp B r /A , K) at the cylinders need not vanish. In princi- 
ple, the magnetic couplin g of the fluid to the cylinders might modify the growth of the 
instability (jLiu et al.l 120071 ) . 



Fortunately, at the frequencies relevant to Princeton MRI experiment (100% run, Q = 
^f^l^2 ~ 150 Hz) the skin depth of stainless steel 5 W = A/2/ostcci^Vio ~ 9 cm (cxsteei is 
the conductivity of stainless steel, see the value in Table. [1]), which is much larger than 
the thickness of the steel vessel surrounding the gallium in the experiment, d w ~ 1.0 cm, 
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so that the magnetic field diffuses rather easily into the boundary. If one considers axial 
currents, the gallium and the steel wall act as resistors in parallel; taking into account their 
conductivities and radial thicknesses, one finds that the resistances of steel walls are much 
larger than the resistance of the liquid gallium [Rj : R n : R IH = 65 : 1 : 21; see Fig. [TJfor 
the subscripts]. Therefore, the currents carried by the steel walls can be neglected for the 
toroidal field, so that the B r B v stress at the boundary is expected to be unimportant and 
an insulating boundary condition suffices. 

The linear growth rate and saturated final state based on partially conducting bound- 
ary condition differs slightly from the results based on perfectly insulating boundary con- 
dition with partially conducting walls, which verifies our argument. Hereafter we ignore 
the conductivity of the stainless steel walls, but regard them as insulating materials, that 
is: r\i = T]u = oo. Please note that in this whole paper the rotation profile used is 
[fti, Q 2 , ^3, ^4] = [4000, 533, 1820, 650] rpm (100% run, Q A > Q 2 , Q u Q 2 , ft 3 , and Q A are the 
rotation speed of the inner cylincer, outer cylinder, inner ring and outer ring, respectively) 
(Tabl eHI), under which the simulations predict that the bulk flow is almost ideal Couette 
state (iKageyama et al.ll2004l ). And this is the desired basic state for MRI to grow. However, 
the experimental scans find that only if we choose [f^, Q 2 , ^3, ^4] — [400,53,146,40] rpm 
(10% run, f2 4 < fl 2 ), the bulk flow profile is found to be very close to the ide al Couette stat e 
with the devia tion (the boundary layer) highly confined near both endcaps (IJi et al.l 120061 ). 
Ji et al.l (120061 ) also shows that non-magnetic quasi-keplerian flows at Reynolds numbers up 
to millions are essentially steady, if the boundary effects have been properly considered. The 
Stewartson layer at the junction of t he rings ma y be smoothed by localized circulation and/or 
turbulence from these instabilities (jLiull2008al ). The rotation profile to be used in the mag- 
netized runs would be [Sl u fi 2 , fi 3 , ^4] = [2400,320,876,240] rpm (60% run, Q A < Q 2 ). The 
difference of the rotation profile at the end caps between in the experiment and in the simu- 
lation could possibly be explained by the wobbling of the inner cylinder in the experimental 
runs due to the diffi culties of aligning the inner cylinder perfectly and the gaps between the 
rings and cylinders (IRoach et al.l 120071 ). And these effects are speculated to be more impor- 
tant with higher Reynolds number Re ~ 10 6 — 10 7 as in the experiment. Unfortunately, 
the modern computers and codes can not afford a simulation with Reynolds number as high 
as several million. Although this is not experimentally realizable, the low-Reynolds-number 
results reported in this paper without wobbling, gaps and starting from an initially ideal 
Couette state with fully insulating magnetic boundary conditions still gives us some useful 
hints about what is going on in the experiment. The result maybe is even closer to what 
is happening in the experiment based on the purely hydrodynamic experiment (steady flow, 
mostly ideal Couette state in the bulk), although this argument needs to be proven. We are 
addressing the problem step by step, adding the complexity into the code one by one. To 
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include all subtle points at one time would perplex the problem and make it difficult for us 
to catch the essence of the problem. The detailed study of the influence of the wobbling and 
gaps between the rings is on the way to investigate the discrepancies between the experimen- 
tal res ults and simulations. Also we want to point out: as demonstrated by iGoodman fc Ji 
( 120021 ). the viscosity of liquid metals is so small as to be almost irrelevant to MRI, at least 
in the linear regime. Thus a low-Reynolds-number run would be fine, at least in this sense. 

Most of results presented in this paper are the results of 100% run except stated explic- 
itly. The simulation results have proved very valuable in design, operation and understanding 
of the experiment. The outline of this paper is as follows. We present the results of the linear 
MRI in Sec. |2j And the results of the nonlinear saturation of MRI are discussed in Sec. |3j 
Sec. |H summarizes the results and presents the final conclusions. 



2. Linear Phase 

2.1. Linear Growth Rate Reduced by the Residual Magnetized Ekman 

Circulation 

The first convincing evidence of the existence of MRI is its linear growth rate (Fig. H|). 
We find this linear growth rate is reduced from 33.1 s _1 with ideal Couette state at both 
endcaps (left panel of Fig. HI the end effect is removed completely by enforcing ideal Couette 
state at both endcaps, thus no Ekman circulation is present) to 21.7 s _1 with two rings 
just as in the real experiment (right panel of Fig. H|). This is due to the resid ual magneti c 



Ekman circulation, which modifies the background flow. Also as pointed out in lLiul (l2008bl ). 
the linear growth rate of an absolute instability in a bounded geometry is reduced by the 
"absorping" end plate, with the reduction extent proportional to 0(T~ 2 ), where the aspect 
ratio is T = h/(r 2 — rx) and h is the height of the cylinders (see Fig. [T]). 

Interestingly the linear growth rate with ideal Couette state at both endcaps (right 
panel of Fig. HJ) matches the growth rate of the faste st growing mode, whi ch naturally is the 



dominant mode, got from a linear code developed by (IGoodman fc Jill2002l ). in which vertical 
periodicity is assumed, but the radial equations are solved directly by finite differences with 
perfectly insulating boundary conditions (Table [2]). The agreement suggests that viscous 
effects are slight, since the no-slip conditions on the horizontal velocities at the endcaps 
differ from the periodicity imposed by the linear code. 
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2.2. Strong Magnetic Field Suppresses MRI with Two Split Rings 

MRI essentially is a weak field instability. It is characteristic of MRI for a strong 
magnetic field to suppress this instability. Our experimental facility can only allow I v < 
1200 A. Thus we need to try carefully to find one set of proper parameters under which one 
growing mode is present with low I v while absent with larger I v . Fig. demonstrates this 
property for simulations in which the boundary rotation rates (Qi, f^, ^3, ^4) are scaled to 
45% of their designed values (TableJT]). 



Nonlinear Saturation 



3.1. Inflowing "Jet" Observed near the Middle Plane 

For Re = 6400, the final state is not steady. Typical time averaged flow and field 
patterns are shown in Fig. The poloidal flux and stream functions are defined so that 



V P = V r e r + V z e z 



r 'e B xV$, 



Bi 



r 'e B xV$, 



which imply V • Vp = and V • Bp = 0. The divergence of the velocity field in these 
compressible but subsonic simulations is nonzero but small. 



We note that the induced toroidal field is around 6% of the initial axially imposed 
magnetic field at this magnetic Reynolds number Re m : -B^max ~ 0.065oo. The most striking 
feature is the inflowing "jet" centered near z = 13.95 cm in Fig. 6] (see also Fig. [7j), which 



is opposite to the usual Ek man circulation (Kageyama et al.l 120041 ). It seems that the rapid 
outflowing "jet" found in (ILiu et al.l l2006bl ) with vertically periodic boundary condition, 
where the position of the "jet" is arbitrary, is shifted to the boundary layer, near the end 
caps in the experiment. 

The inwar d "je t " near the inner cylinder part is a direct consequence of MRI. For MRI- 
stable regimes, iLiul (j2008al ) discuss the influence of independently controlled split endcaps 
upon the poloidal circulation and the influence of the axial magnetic field and Reynolds 
number upon the Stewartson and Ekman layers. In such regimes, an outflowing jet occurs 
near the inner cylinder, in our apparatus, as it does in purely hydrodynamic Ekman circu- 
lation (i.e., A = 0). However in the MRI- unstable regime for the same Elsasser number A, 
the middle clock-wise cell enlarges both horizontally and vertically and eventually dominates 
the other cells, which retreat to the corners if the instability is strong enough. This leads to 
one inflowing "jet" near the inner cylinder (Fig. [7]). 



Further evidence that MRI causes the inflowing jet was found in simulations with the 
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ideal Couette state imposed at the end caps, to remove Ekman circulation (Fig. [8]). Fig. [8] 
shows a clear inflowing jet near the inner cylinder in this case (see also Fig. [9]). This suggests 
that the poloidal circulation seen in the final state in the split-endcap cases (Fig. [H] (b)) is 
caused solely by saturation of MRI rather than (magnetically modified) Ekman circulation. 



3.2. MRI Enhances Radially Outward Angular Momentum Transport at 

Saturation 

Astrophysicists are interested in the angular momentum transport due to MRI since 
MRI is supposed to be the most probable mechanism to explain the fast accretion in the 
astrophysical disks. Fig. [10] displays the time-averaged r-profiles of the radial advective, 
viscous, and magnetic torques, i.e. the angular momentum fluxes integrated over cylinders 
coaxial with the boundaries: 

r a dvcctive,r(r) = 2n dz pr 2 v r v 9 , (2) 







h 



r m agnotic,r( r ) — 2lT j dz ( ^ - ) , (3) 





h 

Fyiscous,r('") 27T / dz 







3 ^ /Up 

or \ r 



(4) 

['total. rv) r advective , r \J) ~\~ F m agnctic,r \J) ~\~ F v iscous,r \J) j (5) 

where dr and dz are the radial and vertical cell sizes respectively. 



In contrast to the final state for vertically periodic boundary conditions (ILiu et al. 



2006bj), the total radial torque is not constant with radius. Since our numerical scheme 
conserves angular momentum exactly, we can infer a vertical flux arising from exchange 
of angular momentum with the endcaps (see also Fig. [TT1) . In the statistical steady state 
(nonlinearly saturated final state) the sum of the angular momentum flux on all boundaries 
should be zero, at least on time average. In this case it is fluctuating since the final state is 
not steady given high Reynolds number. From the gradients of the radial torque, we identify 
four Ekman circulation cells: where dT tota \ tT /dr > (< 0) , the fluid is losing (gaining) 
angular momentum at the endcaps and the boundar y-layer flow is therefore radially inward 



(outward). This is consistent with the discussion in iLiul (j2008al ) of the poloidal circulation 
driven by two split rings. The radial magnetic and advective torques vanish at r\ and T2 
because of the boundary conditions but are important at intermediate radii, especially the 
advective. All components of the radial torque are positive, which means that the angular 
momentum is transported radially outwards. 
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Fig. [TT] displays the time-averaged r-distributions of the vertical total angular momen- 
tum flux at both endcaps (z = Q,h): 

pr 

drpr 2 v z v v , (6) 
r magne tic,z(r, z = 0, h) = 2vr / dr ( - T , (7) 

2 9 ( \ 



0,h) 


IT 

= 2tt / 




J ri 


0,h) 


IT 

= 2tt / 

Jri 


0,h) 


IT 

= 2tt / 

J ri 



r v iscous, z (?", z = 0,h) = 2n I dr 



Since at both endcaps, both the advective and magnetic angular momentum fluxes are zero 
due to the boundary conditions, the total angular momentum at both endcaps r to tai,z(^, z = 
0, h) = r advectiveiZ (r, z = 0, h) + r magne t iCiZ (r, z = 0, h) + r viscouSiZ (r, z = 0, h) are simply the 
viscous angular momentum flux r viscouSjZ (r, z — 0, h). 

From these two figures, we can see that the global angular momentum is entering from 
the inner cylinder then most of it is flowing out from the outer cylinder while the rest of it 
is flowing out at the two end caps. Angular momentum is transported radially outward. 

It is very interesting to derive the dependence of the total torque at the inner cylinder, 
which is mainly responsible for driving the rotation, and of the sum of the vertical torques 
at both endcaps, on Reynolds number (Re). From Fig. [12] (a), we infer the following scalings 
(100% run, Re m = 16): 

(1) In the absence of magnetic field (100% run, 100 < Re < 25600), 

r init iai,r(ri) w 2.69 x 10 4 i?e- a691 ; (9) 

These hydrodynamic states (rather than ideal Couette flow) will be used as the initial con- 
ditions for the MRI experiments. The maximum relative error ratio, defined as the ratio of 
the absolute difference of the value from the simulation and the fitted value (Eq. [9]) over the 
former one, is < 0.4%. 

(2) For 1^ = 1000 A and A = 0.36 (100% run, 100 < Re < 25600). 

rnnaL^rO^i-gsxio 4 ^- - 639 . (10) 

This corresponds to the final states of the MRI experiment. The maximum relative error 
ratio is < 0.3%. From these two scaling laws, the MRI indeed enhances the angular momen- 
tum transport at saturation, though slightly (see discussion below). The enhancement is far 
beyond the arrange of the error bar (see the error bars on panel (a) of Fig. [T2|) . It is reason- 
able for r r in both magnetized and unmagnetized cases to decrease with increasing Re since 
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the viscous coupling to the walls scales as -Re -1 / 2 in standard unmagnetized Ekman layers 
(i.e. the endcaps are corotating solidly with the outer cylinder). In this more complicated 
case (with split rings), we also expect a similar relationshi p, though the magnetized Ekman 



and Stewartson layers complicate the problem (lLiull2008af ) 



If Eq. [9] and Eq. [10] also work at larger Reynolds numbers, so that they may safely be 
extrapolated to the experimental Reynolds number (Re ~ 1.15 x 10 7 , 100% run), then the 
total radial torque of the initial and final state at the inner cylinder may be as large as 

r ini tiai,r(ri) « 0.359 [Newton m] (11) 

and 

r fina i jr (n) w 0.611 [Newton m] (12) 

respectively. Thus, the ratio of the increase of the torque over the initial torque is: (0.611 — 
0.356)/0.356 = 72%, which is quite measurable and indicates that at the experimental 
Reynolds number MRI would dominate the residual magnetic Ekman circulation in the point 
of view of transporting the angular momentum. There are, however, reasons for caution in 
accepting this estimate. For example, the experimental flow may be three-dimensional and 
turbulent, which might result in an even higher torque in the final state, and the absolute 
values of both the exponents seem to decrease at larger Reynolds number and the difference of 
these two exponents is small. These concerns all make the extrapolation of Eq. and Eq. [TO] 
to the experimental Reynolds number a bit risky. Nevertheless, we expect a noticeable torque 
enhancement in the MRI-unstable regime. 

From Fig. [12] (b), we can see that: (1) At larger Reynolds number, (T z (z = 0) + T z (z = 
h))/T r (r = ri) is increased, which means that larger part of the total angular momentum 
is transported vertically. (2) In the MRI stable regime (Re < 1600), the magnetic field 
enhances the vertical transport of the angular momentum. This is also reasonable since the 
magnetic field would align the flow, thus having the cells elongating and penetrating deeper 
into the bulk. The size of the middle cells is increased vertically by the r esid u al mag netic 



Ekman circulation. These are consistent with the conclusions deduced in iLiul ((2008a]). (3) 
In the MRI unstable regime (Re > 3200), the onset of the MRI results in more angular 
momentum transported radially outwards and less vertically. The MRI would increase the 
scale of the middle cell horizontally. Therefore it transports more angular momentum radially 
outwards. 



-10- 



4. Conclusions 



In conclusion of this paper we have simulated the nonlinear development of magnetorota- 
tional instability in a nonideal magnetohydrodynamic Taylor-Couette flow. The simulations 
mimic an on-going experiment except that the conductivity of the stainless steel walls is 
neglected and the simulation is started from an ideal Couette state rather than an actual 
hydrodynamic statistical steady state driven by split end caps. We have also restricted our 
study to smaller fluid Reynolds number (Re) than in the experiment, however we have used 
exactly the same magnetic Reynolds number (Re m ). MRI grows from small amplitudes at 
rates in good agreement with linear analyses without the end cap effects. 

Concerning the MRI simulations with two split independently rotating rings like the 
real experimental facility, we draw the following conclusions: 



1. A clear linear phase is observed; the linear MRI growth rate is reduced by the residual 
magnetized Ekman circulation. 

2. Strong magnetic field suppresses MRI. 

3. In the final state one inflowing; "jet" opposite to the usual Ekman circulation "jet" 
(IKageyama et al.ll2004l ) is found near the inner cylinder, a direct consequence of MRI 
rather than the residual Magnetic Ekman circulation (100% run). 

4. The MRI enhances the angular momentum transport at saturation. (100% run). 

5. The final state contains horizontal fields about 6% as large as the initial vertical field 
for Re m « 20 (100% run). 



We emphasize that these conclusions are based on axisymmetric simulations restricted 
to the range 10 2 < Re, Re m < 10 4 4 with the idealizations mentioned above. The simulation 
results considering the conductivity of the steel container, starting from an actual hydrody- 
namical equilibrium and the comparison with the experimental results would be given in a 
future paper. 
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James Stone for the advice on the ZEUS code and thank Stephen Jardin for the advice 
of implementing full insulating boundary condition. This work was supported by the US 
Department of Energy, NASA under grants ATP03-0084-0106 and APRA04-0000-0152 and 
also by the National Science Foundation under grant AST-0205903. 



- 11 - 



REFERENCES 

Balbus, S. & Hawley, J. 1998, Rev. Mod. Phys., 70, 1 

Burin, M., Schartman, E., Ji, H., Cutler, R., Heitzenroeder, P., Liu, W., Morris, L., & 
Raftopolous, S. 2006, Exp. Fluids, 40, 962 

Fleming, T. P., Stone, J. M., & Hawley, J. F. 2000, Astrophys. J., 457, 355 

Gilman, P. 1971, Phys. Fluids, 14, 7 

Goodman, J. & Ji, H. 2002, J. Fluid Mech., 462, 365 

Hayes, J. G, Norman, M. L., Fiedler, R. A., Bordner, J. O., Li, P. S., Clark, S. E., ud Doula, 
A., & Low., M.-M. M. 2006, Astrophys, J. Suppl, 165, 188 

Hollerbach, R. & Fournier, A. 2004, in MHD Couette Flows: Experiments and Models, 
American Inst. Phys. Conf. Proc., ed. G. Rosner, G. Rudiger, & A. Bonanno, Vol. 
733, 114-121 

Hollerbach, R. & Rudiger, G. 2005, Phys. Rev. Lett., 95, 124501 

Ji, H., Burin, M., Schartman, E., & Goodman, J. 2006, Nature, 444, 343 

Ji, H., Goodman, J., & Kageyama, A. 2001, Mon. Not. R. Astron. Soc, 325, LI 

Kageyama, A., Ji, H., Goodman, J., Chen, F., & Shoshan, E. 2004, J. Phys. Soc. Japan., 
73, 2424 

Knobloch, E. k Julien, K. 2005, Mon. Not. R. Astron. Soc, 000, 1 
Liu, W. 2008a, accepted by PRE 
— . 2008b, in preparation 

Liu, W., Goodman, J., Herron, I., & Ji, H. 2006a, Phys. Rev. E, 74, 056302 
Liu, W., Goodman, J., & Ji, H. 2006b, Astrophys. J., 643, 306 
-. 2007, Phys. Rev. E., 76, 016310 

Noguchi, K., Pariev, V. I., Colgate, S. A., Beckley, H. F., & Nordhaus, J. 2002, Astrophys. 
J., 575, 1151 

Priede, J., Grants, I., & Gunter, G. 2007, Phys. Rev. E, 75, 047303 



- 12 - 



Roach, A., Ji, H., Liu, W., & Goodman, J. 2007, in APS-DPP07, abstract BP8.00084 
Rudiger, G. & Hollerbach, R. 2007, Phys. Rev. E, 76, 068301 

Riidiger, G., Hollerbach, R., Schultz, M., & Shalybkov, D. 2005, Astron. Nachr., 326, 409 

Rudiger, G., Hollerbach, R., Stefani, F., Gundrum, T., Gerbeth, G., & Rosner, R. 2006, 
Astrophys. J., 649, L145 

Sano, T. & Inutsuka, S. 2001, Astrophys. J., 561, L179 

Sisan, D. R., Mujica, N., Tillotson, W. A., Huang, Y., Dorland, W., Hassam, A. B., Anton- 
sen, T. M., & Lathrop, D. P. 2004, Phys. Rev. Lett., 93, 114502 

Stefani, F., Gundrum, T., Gerbeth, G., Rudiger, G., Schultz, M., Szklarski, J., & Hollerbach, 
R. 2006, Phys. Rev. Lett., 97, 184502 

Stefani, F., Gundrum, T., Gerbeth, G., Rudiger, G., Szklarski, J., & Hollerbach, R. 2007, 
New J. Phys., 9, 295 

Szklarski, J. 2007, Astron. Nachr., 328, 499 

Umurhan, O. M., Menou, K., & Regev, O. 2007a, Phys. Rev. E, 76, 036310 
-. 2007b, Phys. Rev. Lett., 98, 034501 

Velikhov, E. P., Ivanov, A. A., Lakhin, V. P., & Serebrennikov, K. S. 2006, Physics Letters 
A., 356, 357 



This preprint was prepared with the AAS IATgX macros v5.2. 



13 



,dw, 



IV 



dw 



®I 



9 



IV 



in in 



ii 



IV 



8.01 cm 
4.32 cm 



Q., 



£2? 



6.165 cm 



R 



IV 



III III 



II 



IV 



4.11 cm 



r l 



IV 



r 2 

29.16 cm 



Fig. 1. — Computational domain for simulations of Princeton MRI experiment. Region (I): 
Inner steel cylinder, angular velocity Q±, magnetic resistivity rjj. (II): outer steel cylinder, 
^2; Vn- liquid gallium, i]Q a ; (IV): vacuum. Thick dash line: insulating inner ring, f2 3 . 

Thick dash-dot line: insulating outer ring, f2 4 . Thin dash line: middle plane. Dimensions: 
radius of the inner cylinder, r\] radius of the outer cylinder, r 2 ; height of the cylinders, h; 
thickness of inner and outer steel cylinder, d w = 0.9525 cm; black rectangles, external coils. 
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Fig. 2. — Figure Courtesy of Mark Nornberg. Experimental wave form of the coil current. 
Left panel: 1^ = 1000 A; right: I v = 400 A. Shown here is the total current in the general 
circuit. Since the six coils are split into two sets of three in parallel, with the upper three in 
series and the bottom three also in series, the current in the coils should be the one in the 
general circuit divided by 2. 




Fig. 3. — Simulated Current pattern in the external coils; the ramp time t coi i = 0.2 s; except 
where stated explicitly, the steady current I v = 1000 A. 
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Fig. A.— 100% run (MRI unstable) and I v = 1000 A. B r . vs. time for Re = 6400, Re m = 20 
sampled outside the fluid at z — 13.95 cm, r = 25.0 cm. "Bottom end cap" is located 
at z — 0. Height h = 27.9 cm. Left panel: endcaps split into two rings, Growth rate 
7 = 21.7 s -1 ; right panel: ideal Couette state at both endcaps, Growth rate 7 = 33.1 s -1 . 



Dimensions 


r± — 7.1 cm 


r 2 = 20.3 cm 


h = 27.9 cm 


d w = 0.9525 cm 


Material Properties 


p Ga = 6.35 g cm" 3 


r] Ga = 2.43 x 10 3 cm 2 s" 1 


Tysteel = 5.73 x 10 3 cm 2 s- 1 


asteei = 1.25 x 10 16 s- 1 


Rotation Profi 


e (100% run) 


1V27T = 4000 rpm 


Q2 /2tt = 533 rpm 


fi 3 /27r = 1820 rpm 


^4/2^ = 650 rpm 


Rotation Profi 


le (45% run) 


Qi/2n = 1600 rpm 


tt 2 /2n = 239.85 rpm 


tt 3 /2ir = 819 rpm 


Ot/27r = 292.5 rpm 



Table 1: Parameters in Gaussian units used in the simulations 
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Fig. 5. — 45% run. B r . vs. time for Re = 6400, Re m = 7.3 sampled outside at z — 13.95 cm, 
r = 25.0 cm. "Bottom end cap" is located at z — 0. Height h = 27.9 cm. left panel: 
lip = 750 A, growth rate 7 = 1.3 s -1 ; right panel: I v = 1200 A, stable. 



Re m 


Re 


n 


Prediction [s 


Simulation [s x ] 


16 
16 


6400 
6400 


1 
2 


33.7 
13.8 


33.1 



Table 2: 100% run. Growth rates from semianalytic linear analysis vs. simulation. 
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Fig. 6. — 100% run. Contour plots of final-state velocities and fields (MRI unstable). Re = 
6400, Re m = 20. I v = 1000 A, A = 0.36. (a) Poloidal flux function * (Gauss cm 2 ) (b) 
Poloidal stream function <3> (cm 2 s _1 ) (c) toroidal field B v (Gauss) (d) angular velocity Q = 
r _1 V^ (rads" 1 ). In panel (b), dashed line indicates the clock-wise poloidal circulation while 
solid line indicates the anti-clock-wise poloidal circulation. 
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Fig. 7. — Corresponding to Fig. Efb). 100% run (MRI unstable). Time averaged v r vs. 
radius r on the midplane (z = h/2). Re = 6400, Re m = 20. I v = 1000 A, A = 0.36. 
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Fig. 8. — 100% rotation run (MRI unstable). Contour plots of poloidal stream function 
$ (cm 2 s _1 ) for Re = 6400 with ideal Couette state at the end caps. 1^ = 1000 A, A = 0.36. 
Dashed line indicates the clock-wise poloidal circulation while solid line indicates the anti- 
clock-wise poloidal circulation. 
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Fig. 9. — Corresponding to Fig. El 100% rotation run (MRI unstable). Time averaged v r vs. 
radius r on the midplane (z = h/2) for Re = 6400 with ideal Couette state at the end caps. 
I v = 1000 A, A = 0.36. 
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Fig. 10. — 100% run (MRI unstable). Time averaged ^-integrated radial angular momentum 
fluxes versus r at saturation. Re = 6400, Re m = 20. 1^ = 1000 A, A = 0.36. Dash line , 
viscous torque; dash dot line , magnetic torque; long dash line , advective torque; solid line, 
total torque. 
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Fig. 11. — 100% run (MRI unstable). Time averaged vertical total angular momentum flux 
at both endcaps versus r at saturation. Re = 6400, Re m = 20. I v = 1000 A, A = 0.36. 
Solid line: bottom endcap (z — 0); dashed line: top endcap (z — h). 
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Fig. 12. — (Color) (a) Total radial torque at the inner cylinder (b) Sum of the total vertical 
torques at both endcaps versus Re. Note that in the simulations the magnetic diffusivity rj is 
fixed to the experimental value rj = 2430 cm 2 s _1 (TableJIJ), however the kinetic viscosity is 
varied for the purpose of extrapolation. In both panel, red colors: MRI unstable; black color: 
MRI stable. In panel (a), dashed lines have slopes of —0.691 (initial state) and dash-dot 
lines —0.639 (final state). The error bars on panel (a) (very small, typically ~ 1%, hard to 
see on panel (a)) are the standard time-deviations. If Re < 3200, the error bars are zero 
in both magnetized and unmagnetized cases, which suggests the system is steady with low 
Reynolds number. 



